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Abstract. The evolution of optically selected quasars is 
usually supposed to be well described by a single con- 
stant evolution parameter, either or kjj, depending 
whether we refer to luminosity or density evolution. In 
this paper we present a study of the variations of the 
evolution parameters with rcdshift, for different cosmo- 
logical models, in order to probe the differential evolution 
with redshift. Two different quasar samples have been an- 
alyzed, the AAT Boyle's et al. and the LBQS catalogues. 
Basically, these samples are divided in redshift intervals 
and in each of them k^ and ko are estimated by forc- 
ing that (V/V max ) — 0.5. The dependence with respect 
to the cosmological parameters is small. Both AAT and 
LBQS show roughly the same tendencies. LBQS, however, 
shows strong fluctuations, whose origin is not statistical 
but rather due to the selection criteria. A discussion on se- 
lection techniques, biases and binning effects explains the 
differences between these results. We finally conclude that 
the evolution parameter is almost constant in the redshift 
range 0.7 < z < 1.7, at least within 2a, while it decreases 
slightly afterwards. Results depend on the binning chosen 
(but not in a very significant way). The method has been 
tested with Monte-Carlo simulated catalogues in order to 
give a better understanding of the results coming from 
the real catalogues. A correlation between k^ (&£>) and 
(V/V max ) is also derived and is used for the calculation of 
the error bars on the evolution parameter. 

Key words: Cosmology: observational tests - Quasars: 
general, evolution 



1. Introduction 

Since the first application of the V/V m ax test by Schmidt 
in 1968, the luminosity function of quasars is known to 
undergo a strong evolution, in the sense that the density 
of the most luminous quasars was far higher in the past. 
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This was firstly interpreted in terms of either Pure Den- 
sity (PDE - Schmidt, 1968) or Pure Luminosity Evolution 
(PLE - Mathez, 1976), and phenomenologically modeled 
with a single free parameter law. More complex models 
began to appear a few years later, such as a Luminos- 
ity Dependent Density Evolution (LDDE), proposed by 
Schmidt and Green in 1983. However, there wasn't any 
privileged model until the late 80's, when Boyle et al. 
(1988) favored PLE as a best fit of the observed lumi- 
nosity function of a UVX sample. A similar recent result 
is found for X-Ray selected samples of quasars (Page et 
al. 1997; Jones et al., 1996; Boyle et al., 1993). The cur- 
rently observed evolution of quasars appears however to 
be more complex nowadays, since the first indication of 
a reversing evolutionary trend around a redshift in the 
range [2.5,3] was observed (Shaver, 1994; Schmidt, Schnei- 
der and Gunn 1995; Warren, Hewett and Osmer 1994; Pei, 

1995) . A steepening of the luminosity function towards 
high redshifts is advocated by Goldschmidt and Miller 
(1997). 

Meanwhile, the evolution of the cosmic star formation 
density was found to be strikingly similar to that shown 
by Shaver for QSOs, although the maximum of the star 
formation density is attained at a lower redshift (Madau, 

1996) . In spite of this difference, the similarity of the vari- 
ations of QSO luminosity density and of the field galaxy 
star formation has been interpreted as a clue that both 
phenomena could be closely linked (Boyle and Terlevitch, 
1998; Silk and Rees 1998). Furthermore, the similarity is 
even more striking (both curves have a maximum around 
the redshift z=2.5) after applying the necessary correction 
for dust extinction in the density of high redshift galaxies 
(Shaver et al., 1998). This correction may be as high as 
a factor of ten since these galaxies are observed in their 
rest-frame ultraviolet. 

In the same time, complex models of this evolution 
begin to appear: 

• At high redshift, the growth of Dark Matter ha- 
los according to Press-Schechter formalism and the paral- 
lel Eddington-limited growth of accreting Massive Black 
Holes (MBH) (Haiman and Loeb, 1997) are likely to in- 
duce a decrease of density with increasing redshift ('neg- 
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ative DE') (Haiman and Loeb, 1997; Cavaliere and Vit- 
torini, 1998; Krivitsky and Kontorovich, 1998; Novosyad- 
lyj and Chornij, 1997). 

• Around a redshift z=2.5, there could be two phenom- 
ena: a transition from high to low efficiency in advection- 
dominatcd flows, followed by the decline of accretion rate, 
giving luminosity proportional to (1 + z) k with k slightly 
variable around 3 (Yi, 1996). Alternatively, growing galax- 
ies could assemble in groups of typically 5 10 12 Mq, where 
tidal effects refuel the MBH at lower and lower rate, trans- 
lating into 'positive LE' (Cavaliere and Vittorini 1998). 

• At intermediate redshift, viscous instabilities in- 
duce long term, high amplitude variations of the ac- 
cretion rate. The fraction of time spent at each 
luminosity level, convolved with the mass distribu- 
tion gives the luminosity function and its evolution 
(Siemiginowska and Elvis, 1997). Galaxy collisions pro- 
vide a mechanism which fuels galactic nuclei with gas in 
dense environments, giving raise to quasars in low lumi- 
nosity galaxies (Lake, Katz and Moore, 1997). 

• In addition, gravitational lensing mimics luminos- 
ity evolution in flux-limited samples, amplifying the more 
distant quasars, but whose detailed effect is not known, in 
particular the induced bias selection effects in magnitude 
limited quasar samples. 

Van Waebeke et al. (1996) define a new test, the 
V/Vmax statistics, and apply it to the AAT quasar sam- 
ple (Boyle et al., 1990). They show that, contrary to the 
usual (V/V m ax) test, the V/V max statistics leads to con- 
straints on cosmological parameters. However these results 
rely on an arbitrary model of PLE, since the evolutionary 
effects dominate the cosmological effects in the QSO dis- 
tribution, and we have to be sure of the reliability of the 
evolution model before constraining the cosmological pa- 
rameters. So, the analysis in Van Waebeke et al. (1996) is 
rather a test of compatibility between an evolution model 
and a couple of cosmological parameters and there is a 
need for an independent better understanding of the QSO 
evolution before applying such cosmological tests. 

A number of projects aim at assembling complete 
and homogeneous samples of several thousands of quasars 
which will allow more reliable and more subtle analyses, 
leading, thus, to verifications or improvement of complex 
theories on quasar evolution. 

Finally, there is some hope for progress in the under- 
standing of the evolution of the stellar formation rate and 
of quasars, and in the determination of cosmological pa- 
rameters compatible with the observed quasar distribu- 
tion. Note that this is a necessary framework in order to 
make a precise estimate of the effects of the re-ionization 
in the future temperature maps of Planck-Surveyor. 

The aim of this paper is to explore in details the pos- 
sibility for the evolution to depend on the redshift, and 
on whether we can distinguish between artefact and phys- 
ical evolution. For the present study we decided to adopt 
a cosmology and an evolution model with a single evolu- 



tion parameter which is allowed to depend on the redshift 
(and/or the magnitude in Section 4). This should lead to 
more complex evolution laws, which may help in theoret- 
ical understanding of quasar evolution and in turn could 
be the basis of future, more realistic cosmological tests. 

The technical basis of this work is the (V/V ma x} test, 
performed in bins of redshift. In §2 and 3 we make a brief 
description of the quasar samples and the evolution mod- 
els used. In §4 we give a justification of our choice to study 
quasar evolution in redshift bins and in §5 we present the 
different redshift binning modes used for our study. §6 con- 
tains the results obtained from the simulated catalogues. 
§7 contains the results for the real quasar samples. More 
precisely, it shows the measured dependence of the evolu- 
tion parameter versus redshift, for a few sets of cosmolog- 
ical parameters. Finally, in §8 we discuss the main results 
of our study. 

2. Samples used 

In this study two complete quasar samples have 
been used, the AAT catalogue of Boyle et al. 
(Boyle et al., 1990) and the Large Bright Quasar Sur- 
vey catalogue (Hewett et al., 1995), hereafter AAT and 
LBQS respectively. The AAT catalogue contains 400 faint 
quasars selected by their UV excess, it is claimed to be 
complete within 18 < to < 21 and 0.2 < z < 2.29. LBQS 
consists of 1055 optically selected quasars of intermediate 
brightness, with 16.5 < vhbj < 19, within the redshift 
range 0.2 < z < 3.4. Various methods have been used for 
the selection of this catalogue, combining a color selection 
technique and slitless spectroscopic data. LBQS is almost 
three times bigger in size and much more profound in red- 
shift than AAT, but its selection criteria depend on the 
redshift. This means that the selection biases differ from 
one bin of redshift to another, which will lead to dramatic 
effects on the measured evolution. In the case of AAT, 
the UVX selection technique applied in the whole redshift 
range insures the fact that the biases can be comparable 
throughout the whole sample. The differences between the 
two samples may influence the results on quasar evolution 
thus a direct comparison may not be possible. 

3. Evolution models and evolution parameter 

3.1. PLE and PDE 

In order to explain the observed non-uniform spatial dis- 
tribution of quasars, several empirical evolution models 
have been proposed such as the PLE (Pure Luminos- 
ity Evolution) model, the PDE (Pure Density Evolution 
model) and the LDDE (Luminosity Dependent Density 
Evolution). In this paper, we will make an extensive use 
of PLE and PDE, whose introduction was purely phe- 
nomenological. PLE makes the hypothesis of a constant 
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object space density and supposes the following form for 
the luminosity: 



L(z) = L(0)e(z) 



(1) 



where e(z) is an one parameter evolution law. On the con- 
trary, PDE assumes a constant quasar luminosity and a 
redshift-dependent space density: 



P (z) = p(0)e(z) 



3.2. Power Law and Exponential Parametrizations 



(2) 



The most commonly used expressions for e(z) are the 
power low form 



e{z) = (l + z) k 

and the exponential form 

e(z) — exp(kT(zj) 



(3) 



(4) 



where r(z) is the lookback time, and k is the evolution 
parameter (k^ and ko in the cases of a PLE and a PDE, 
respectively) . 

The basis of this work is the (V/V max ) test, or rather 
the (V e /Va) test (Avni and Bahcall, 1980), with all vol- 
umes computed according to Mathez et al., 1996. 

In the following we examine the variations of the 
evolution parameter in magnitude and redshift bins for 
PLE and PDE models. All kj, (&d) calculated here are 
the optimum values which insure a (V/V max ) = 1/2 
{{W/W m ax) = 1/2) 1 within each bin, a necessary (but not 
sufficient) condition for a uniform V/V max (W/W max ) dis- 
tribution within the interval [0, 1]. The dependence of fc^'s 
(ko's) versus redshift is measured and all calculations are 
made for several cosmological models. During our work, 
we examined cosmological models with fio and A varying 
between and 1 with steps of 0.25. In the present paper we 
chose to demonstrate our results for the flat univers cos- 
mological models (Slo + A = 1), excluding the rather non- 
realistic model (fi ,A) = (0,1). In all plots representing 
fc^ or ko versus z, the values of the evolution parameter 
increase with increasing Qq and with decreasing A, unless 
specified otherwise. So the lower curve corresponds to the 
model (fio, A) = (0.25,0.75) while the upper one to the 
model (fio: A) = (1.0,0.0). The above method allows a di- 
rect determination of the evolution parameter error bars; 
the la region of {V/V max ) 1 calculated assuming no evolu- 
tion, is given by 1/2 ± la with a = (vl2iV) _1 , where N is 
the number of quasars within each bin. A corresponding 
error bar on the evolution parameter can thus be derived 
from this la interval of {V/V max ). We will go back to this 
point in the Appendix. 

Our method differs significantly from that used by 
Boyle et al., 1988 on the AAT sample. They showed that 



PLE power law with kh — 3.5 is a good fit to the AAT 
sample, with a technic which consists in calculating the 
luminosity function within redshift bins and in fitting the 
results to a luminosity function model. This method has 
the drawback to consider the evolution parameter as a 
free adjustable parameter, as well as the slopes of the lu- 
minosity function, the density, and the absolute magni- 
tude linking the two slopes. It docs not take into account 
the physical information that the quasars should be uni- 
formly distributed, once the evolution has been corrected. 
Indeed, this information allows the evolution parameter to 
be determined directly from a V/V max test without any 
hypothesis on the luminosity function. 



4. Evolution in bins of magnitude 

4-1. Apparent magnitude 

In case of PLE it is very difficult to divide a sample in 
bins of absolute magnitude, even through complex meth- 
ods, because absolute magnitude has to be estimated at 
a given epoch and thus depends on the evolution law. On 
the contrary, a division in bins of apparent magnitude, m, 
is very easy to operate. The AAT sample has been divided 
into 5 bins which contain equal numbers of quasars. In Fig. 
1 we present the evolution parameter ko versus m under 
a power law hypothesis. We notice a diminishing trend up 
to m ~ 20.5 reversing towards fainter magnitudes. This 
trend can be a priori understood in terms of variations of 
evolution with either redshift or absolute magnitude since 
both are more or less correlated with apparent magnitude. 




19.5 20.0 

Apparent maqrkjde 



Fig. 1. Variations of the evolution parameter fcn with appar- 
ent magnitude for fiat universe cosmologies, with a power law 
evolution model. The different curves correspond to the differ- 
ent cosmological models specified in the text. 



W(z) is defined as: dW(z) = p(z)dV(z) 
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4-2. Absolute magnitude 

Analyzing a sample of quasars with a PDE hypothesis, and 
with Mo < —23 for the cosmological models with qo = 0.1 
and go = 0.5, Schmidt and Green (1983) derived that the 
evolution parameter ku increases with absolute magni- 
tude. Making a similar analysis on AAT and LBQS cata- 
logues (which span a larger absolute magnitude range), we 
find no clear tendency for the evolution parameter ku ver- 
sus absolute magnitude, whatever the cosmological model 
used. The samples have both been divided in 5 bins of ab- 
solute magnitude with a constant number of objects (~ 80 
for AAT, - 210 for LBQS). Fig. 2 illustrates the results for 
the zero curvature cosmological models examined above. 
On the right side of the figures we present a typical error 
bar. Note that the scaling of the vertical axis is the same 
in both Figs. 1 and 2. 

From this previous analysis, we conclude that Fig. 1 
could reflect a redshift dependence, not an absolute mag- 
nitude effect, and that is why this paper is devoted to the 
redshift dependence of the evolution. 



the evolution laws derived for quasars, the redshift enters 
through the factor (1 + z) k (sec equation 3). With such 

a law and binning (Bl) in Table 1, the ratio ^ ) 

which governs the evolution rate within bin decreases 
strongly from bin to bin and this is likely to bias the anal- 
ysis. If we want to take into account the dependence of the 
evolution rate with redshift, we must impose the redshift 
intervals' limits, Zi, such as: 

1 + Zj+l _ 1 + Zi+2 

l + Zl ~ 1 + z l+1 

With this binning (B2), the evolution should be compara- 
ble within all bins and we should be expecting a somehow 
more uniform (V/V max } distribution, provided that k does 
not vary too much from bin to bin. 



Table 1. Bins' redshift limits and quasar numbers for AAT 
sample. The Zi for the third binning have been computed for 
the cosmological model (f2o,A) = (0.5,0.5). 



5. Bins of Redshift 

The quasar samples have been divided into several redshift 
bins in order to examine the redshift dependence of the 
evolution parameter and the influence of the cosmological 
model. The AAT analysis is limited to the redshift interval 
[0.3,2.2], because, at low and high redshift, the sample is 
obviously not complete. The binning intervals are defined 
here. 

The binning was chosen so as to be optimized with re- 
spect to different criteria. The first one is that bins should 
contain a sufficient number of quasars in order to make 
a reliable statistical analysis. The second one is that they 
should be narrow enough so that the evolution parameter 
may be considered as constant within a single bin. The 
third criterion is that, according to the parametrization 
of the evolution hypothesis, we must chose the binning in 
order to distribute the amount of evolution equally within 
bins. 

Our first choice of binning insures equal numbers of 
quasars, in order to distribute equally the statistical noise 
between bins and to minimize it, thus satisfying the first 
criterion above. The AAT sample was divided in 5 bins of 
about 80 quasars each (binning (Bl), Table 1) while the 
LBQS sample was divided in 11 bins each one containing 
about 100 quasars. Redshift limits for LBQS are: 0.20, 
0.37, 0.57, 0.75, 0.96, 1.14, 1.36, 1.81, 2.10, 2.52, 3.36. 

We must be sure however that the variations of the 
observed evolution we are looking for are not biased by 
the binning chosen. It could be the case if the evolution 
rate differed too much from bin to bin. This is why, for 
each assumed evolution law e(z), we look for a binning 
ensuring equal ratios e(z i+ i) / e( Zi ) . In the case of a power 
law parametrization for example, used to describe most of 



(Bl) 


Zi 


0.30 


0.86 


1.22 


1.57 


1.87 


2.2 




n 




78 


78 


78 


78 


78 


(B2) 


Zi 


0.30 


0.56 


0.86 


1.23 


1.67 


2.2 




n 




19 


66 


72 


104 


124 


(B3) 


Zi 




0.30 


0.51 


0.80 


1.27 


2.2 




n 






15 


47 


99 


222 



The exponential parametrization leads to another bin- 
ning, also given in Table 1 (binning (B3)) for the cosmolog- 
ical model (f2o,A) = (0.5,0.5). Here, the look-back time 
intervals are constant within bins (note that the bins' lim- 
its depend slightly on the cosmological model). We chose 
to divide the sample in 4 bins because of the very small 
numbers of quasars contained in the first two bins in a 
choice of a division in 5 redshift intervals (13 and 33). We 
remind the reader that a high number of quasars within 
bins is very important because of the statistical nature not 
only of the test, but of the evolution itself, which is a fur- 
ther motivation for large redshift bins. But too large a bin 
may bias the evolution parameter estimation if the evo- 
lution rate changes too strongly with the redshift, which 
may be the case in the last redshift interval for this bin- 
ning. 

6. Tests with simulated catalogues 

The method has first been validated with Monte-Carlo 
simulated catalogues. The aim of these tests was to exam- 
ine the results obtained by analyzing simulated catalogues, 
constructed under specific hypotheses for the whole red- 
shift interval, whose properties are a priori known. Here 
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Fig. 2. LDDE for a power law evolution model, for AAT (left) and LBQS (right): variations of the evolution parameter kn with 
absolute magnitude for 4 flat universe cosmological models: 1) (fi , A)=(0.25,0.75), 2) (fi , A)=(0.5,0.5), 3) (fi , A)=(0.75,0.25) 
and 4) (too, A) = (l. 0,0.0). There is no clear trend for fee with absolute magnitude. The bar on the right gives the error bar's size 
(la. 



we present two different tests applied to simulated cat- 
alogues whose characteristics mimics the properties of 
AAT. The details of the simulation algorithm can be found 
in Mathez et al. 1996. 



For the first test, 38 simulated catalogues of 400 
quasars each, were built under a constant luminosity 
evolution hypothesis, a power law parametrization with 
k L = 3.47, within the redshift region z e [0.3,2.2]. 
The cosmological model used for the construction was 
(fio = 0.5, A = 0.5). Fig. 3 shows the measured evolu- 
tion in redshift bins (with the binning B2) on these cata- 
logues. The solid line indicates the true value fc^ = 3.47, 
the dashed-dotted line includes the 68% of the values and 
the dashed line corresponds to the variations of the me- 
dian value with redshift. The analysis with binning 
(Bl) on the same catalogues gives a similar result, which 
proves that for a constant evolution rate, the method is 
not sensitive to the chosen binning. 



1-2 
<z> 



Fig. 3. kL versus z for 38 simulated catalogues, binning B2. 
The dashed-dotted line includes the 68% of the fez/s values, 
while the dashed line illustrates the variations with redshift 
of the median value. All catalogues were constructed under a 
PLE, power law hypothesis of fez, = 3.47 for the cosmological 
model (to ,A) = (0.5,0.5). 



For the second test, 44 simulated catalogues were build 
under a PLE, exponential parametrization (k^ = 5), for 
the same cosmological model as above. The analysis has 
been made for binning (B3) and assuming a (wrong) PDE 
hypothesis. The results on Fig. 4 shows that a density 
evolution is measured, with k& slightly lower than the 
true k^. It should not be surprising, indeed with a V/V max 
method, it is impossible to separate luminosity and density 
evolution. On the contrary, the distinction between these 
two kinds of evolution is possible with methods based on 
luminosity function hypothesis (as in Boyle et al., 1988), 
but as explained in Section 3, it suffers from some other 
problems. 



7. Evolution parameter versus z 

The evolution parameter is now calculated within each 
redshift bin, by forcing that (V/V max ) = 0.5 (or 
(W/W max ) = 0.5). An error bar is also derived, from the 
ler interval of V/V max , according to the method described 
in details in the Appendix. 

7.1. PLE results 

Fig. 5 shows 1cl(z) for AAT catalogue for power law and 
exponential forms, and zero curvature cosmological mod- 
els. Binning (Bl) is the first to be examined. 
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0.8 1.0 1.2 1.4 1.6 1.8 2.0 

<z> 




0.8 1.0 



Fig. 5. AAT: The sensitivity of fci to cosmological models and its variations with redshift, for equally populated bins (binning 
Bl), for a power law parametrization (left graph) and the exponential parametrization (right graph). The models treated are 
those with zero curvature. 



0.6 0.8 



1.2 1.4 



Fig. 4. kD versus z for 44 simulated catalogues, binning B3. 
The dashed-dotted line includes the 68% of the fco's values, 
while the dashed line illustrates the variations with redshift 
of the median value. All catalogues were constructed under a 
PLE, exponential law hypothesis of /cl = 5 for the cosmological 
model (fio,A) = (0.5,0.5), and analysis is made under PDE 
exponential law hypothesis. 



The left plot corresponds to a power law and the right 
one to the exponential law parametrization. In both cases 
kh decreases slightly with redshift. The same analysis on 
the LBQS (Fig. 6), results in a figure dominated by strong 
fluctuations, absent in the result given by the AAT (Fig. 
5). The fact that AAT and LBQS have not the same mag- 
nitude range is not enough to explain why the evolution 
is so different between these two catalogues. Fig. 6 also 
shows a global decreasing trend of Ul with increasing red- 
shift, which is marginally significant accounting for the 
error bars. 



The reason for this discrepancy should first be searched 
in the incompleteness and/or the selection bias of the cat- 
alogues, in particular in the LBQS which has a complex se- 
lection criteria. The combination used (hyperscnsibilizcd 
IIIa-J+GG395) had a maximum width of 2000 A with 
a maximum response at a wavelength of around 4200 A 
and a blue limit imposed by the atmospheric cut-off at ~ 
3200 A while the red one was due to the emulsion. Ob- 
jects either with a detected emission line or whose median 
wavelength of the SED was found bluewards of the cen- 
tral wavelength were classified as quasar candidates while 
the others were rejected (Hewett et al., 1995). The most 
intense emission lines of quasars in the waveband we are 
interested in along with their relative intensities are listed 
in Table 2. 



Table 2. Quasars' strong emission lines and relative intensities 
taken from current literature. 



Emission 


wavelength (A) 


relative 


line 




intensity 


Lya 


1216 


1 


Civ 


1549 


0.5 


Cm] 


1909 


0.2 


Mgn 


2798 


0.2 



As redshift increases, the quasar strong emission lines 
appear successively in the blue and in the red part of the 
LBQS spectral range, leading each time to a shift of the 
median wavelength. This gives rise first to an excess and 
then to a deficiency of quasar candidates. Fig. 7 gives our 
estimation on this effect on the detection rate, and finally 
in the evolution parameter fcj,. We estimated that the fil- 
ter's response gets its maximum value at a wavelength of 
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Fig. 6. The sensitivity of kh to cosmological models and its variations with redshift for the LBQS analysis in equally populated 
bins (binning Bl). The left graph corresponds to a power law hypothesis; the right one to an exponential form. The models 
treated are those with zero curvature. The fluctuations seen above are most probably due to the passage of strong emission 
lines through the filter's mean wavelength. 




Fig. 7. Relative excess and deficiency of apparent evolution 
due to strongest emission lines in the LBQS. 

~ 3700 A and that it remains almost constant up to a 
A <~ 4900 A. Based on these hypotheses we first calcu- 
lated the redshift regions where each emission line would 
contribute to an excess (and to a deficiency) depending 
on its relative intensity. The conversion of wavelengths to 
redshifts has been made by using the forms: 



(A L + AA L )(1 + z 1L ) = Ai 
(A L ±AA L )(l + z 2L± ) = A 2 
(A L -AA L )(l + z 3L ) = A3 



(6) 
(7) 
(8) 



with Ai = 3700 A , A 2 = 4200 A and A 3 = 4900 A , 
Al the wavelengths given in Table 2 and AAl the Full 
Width Half Maximum (FWHM) of each line. Z\ L is the 
redshift at which the emission line L enters the detection 



range, contributing to an apparent quasar excess. z 2 l+ 
marks the end of this contribution. At z = z 2 l- the emis- 
sion line enters the filter's red region and causes an ap- 
parent deficiency of objects and after z$l it is no longer 
detectable. As FWHM and EW vary from one quasar to 
the other and from one emission line to the other, approxi- 
mate values have been used, taken from Ulrich, 1989. The 
composition of Fig. 7 (whose similarity with Fig. 6 is ob- 
vious), has been made by superimposing the results for 
the 4 emission lines. The fcj,'s maxima approximately co- 
incide with the probable high detection rate points giving 
a possible explanation for the k L 's fluctuations with mean 
redshift. According to Hewett et al. (1995), other factors 
are likely to influence also the quasars' detection rate but 
their analysis is outside the aims of this paper. 

Right afterwards we test binnings (B2) and (B3) for 
AAT sample. Fig. 8 shows the variations of with red- 
shift z for these two cases. The values of kL are almost con- 
stant within error bars in the redshift range 0.7 < z < 1.7. 
However, the values of kL in small redshift bins are not 
very reliable, as we can see from the first error bar in both 
graphs, due to the very small number of quasars (see Table 
1). We repeated the procedure described above but this 
time we divided the sample into 4 bins, so as to have more 
objects (32) in the first bin. We found that the evolution 
parameter was significantly higher in this bin and constant 
in the others but we cannot use these results: because of 
the large size of the three other redshift intervals, the test 
looses its sensitivity. From the above we conclude that a 
PLE (or a PDE as shown later on), with fc^ — 3.5 for a 
power law and kL — 9 for an exponential, fits the sample, 
as already demonstrated by Boyle et al. (1988) (but our 
method differs significantly, as explained in Section 3.2). 
We also applied the analysis to LBQS for the equivalent 
of binnings (B2) and (B3) as well as for random binnings, 
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even though we consider this sample as highly biased. We 
find in average no significant decrease of k with increasing 
z, up to z ~ 2. Fig. 9 illustrates these results for binning 
(B2). Peaks are far less intense but excess and deficiency 
are found in the predicted places (compare with Fig. 7). 




0.5 1.0 1.5 2.0 

<z> 



Fig. 9. fez, variations with redshift for LBQS, for binning (B2). 
fc_L is rather constant within the redshift range 0.3 < z < 2. 



7.2. PDE results 

The same work has been made under the hypothesis of a 
PDE, for both power law and exponential laws. In Figs. 
10 and 12 we present the results for AAT for the three 
binnings and in Fig. 11 the results for LBQS for the first 
binning. 

In all cases we find the same characteristics and the 
same trends as for the PLE hypothesis. As already men- 
tioned in Section 6 about the simulated catalogues, a 
density evolution almost as strong as a luminosity evo- 
lution is measured because of the degeneracy between 
these two kinds of evolution, which is not broken by the 
V/V m ax method. All figures concerning binning (Bl) for 
AAT point out that whatever the hypothesis on the evo- 
lution law, the evolution parameter declines slightly with 
redshift. The similarity between the results of the PLE and 
PDE hypotheses indicates that both evolution models fit 
equally AAT sample. In the case of LBQS the evolution 
parameter's peaks become less intense but they are always 
present. 

Furthermore, the small variations with redshift of kh 
and ku imply that both models, when assuming a con- 
stant evolution parameter, are not far from being correct, 
at least in the redshift range [0.7, 1.7]. With binning (Bl) 
for AAT however, the picture (Fig. 5 and 10) was rather 
of a regular decrease of the evolution parameters with in- 
creasing z. In order to check this point we made the same 
analysis in a series of random binnings, for both AAT and 
LBQS. We found that most binnings rather point to a 



constant evolution parameter, at least up to z = 1.7 and 
to a small tendency towards a regular decrease above this 
redshift. 



7.3. Summary 

We will now give a different representation of our results 
concerning AAT, by calculating the relative luminosity 
and density, normalized at z = 2.2, for both PLE and PDE 
models. Fig. 13 illustrates log L(z)/L{2.2) (solid lines) 
and log p(z) / p(2.2) (dashed lines) versus z. Diamonds and 
squares correspond to binning (Bl) given in Table 1, for 
a power law parametrization. Triangles and stars corre- 
spond to binning (B2), for the same parametrization. The 
two curves without symbols represent the third binning, 
(B3), and an exponential parametrization. Error bars have 
been plotted for each binning. As it may be noticed, they 
become wider as the redshift is smaller because of the 
propagation of errors in the calculation of the normaliza- 
tion constants. All calculations have been made for the 
(flo, A) = (0.5,0.5) model. There is a remarkable similar- 
ity with previous results (see e.g. Fig. 7 of Shaver (1994), 
Fig. 3 of Pei (1995) and Fig. 1 of Shaver et al., 1998, all 
obtained from different procedures) as well as with the 
theoretical results of Yi (1996, Fig. 2). In the inset figure 
we present the same quantities calculated for a constant 
evolution parameter. The solid line represents a power law 
parametrization of k = 3.5 while the dashed line repre- 
sents an exponential parametrization of k = 9. 




Fig. 13. Relative Luminosity (solid lines)/density (dashed 
lines) normalized to z = 2.2 versus z. D,0: (Bl), power law 
parametrization. A,*: (B2), same parametrization. Lines with- 
out plotting symbols: (B3), exponential parametrization. Inset 
figure: the solid line represents a power law parametrization 
of k = 3.5 while the dashed one represents an exponential 
parametrization of k = 9. 
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0.6 0.8 1.0 1.2 1.4 1.6 1.8 0.4 0.6 0.8 1.0 1.2 1.4 1.6 



Fig. 8. fc_L variations with < z > for AAT, for fiat universe models and for binning B2 given by equation 5 (left graph) and for 
binning B3 ensuring equal look-back time intervals (right graph). (Note the difference of the scales of the vertical axis on the 
two graphs.) 



0.8 1.0 1.2 1.4 1.6 1.8 2.0 

<z> 




Fig. 10. ko versus < z > for AAT catalogue, binning Bl. Left graph: power law PDE, right graph: exponential PDE. 




Fig. 11. ko versus < z > for LBQS catalogue, binning Bl. The left (right) graph corresponds to the power law (exponential) 
parametrization. 
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0.4 0.6 



1.4 1.6 



Fig. 12. The analysis of AAT catalogue: fcc versus < z > for flat universe models and for binnings (B2) and (B3) given in 
Table 1, for the flat universe models. 



8. Discussion 

This paper presents an analysis of the evolution rate of 
quasars at various epochs. Our method differs from most 
previous ones in that it docs not make use of any de- 
termination of luminosity function. It solely consists in 
computing the evolution parameters fcz, and fcu (assum- 
ing PLE and PDE respectively) in bins of rcdshift, such 
that the V/V max are uniformly distributed. However, the 
various selection criteria applied in the construction of the 
two catalogues (AAT and LBQS) as well as the different 
rcdshift and magnitude limits make a direct comparison 
of the results quiet difficult. The method has been tested 
with Monte-Carlo simulated catalogues. Three binnings 
have been adopted, the first one with equal numbers per 
bin, the second and the third with a priori similar evolu- 
tion rates inside each bin, according to the evolution law. 
Our results are the following: 

• Both samples roughly show the same large trends, 
however modulations do appear in the LBQS results, 
which are likely to be correlated with the crossing of the 
main emission lines from the blue to the red side of the 
available spectral range. 

• All results on the evolution parameters are quite sim- 
ilar whatever the hypothesis, PLE or PDE, mainly due to 
the inefficiency of the V/V max test to distinguishing be- 
tween density and luminosity evolution. 

• Similar results are also obtained for both power law 
and exponential parametrizations. Certainly a single phe- 
nomenon cannot be described by two different laws but we 
note that there is no significant difference between these 
laws in the redshift range z G [0.7,2.2], as seen in the 
inset plot in Fig. f3. The determination of quasars' evo- 
lution rate at high redshifts will also determine which of 
the proposed models is the appropriate one (if any). 



• In bins with equally distributed evolution, however, 
both parameters k^ and ko show far less variations with 
redshift and within error bars we can suppose that they 
are constant (2.5 < k^ < 4, 1.5 < ko < 3 in a power law 
paramctrization) . 

• fci (ko) and (V/V max ) are linearly correlated, as 
shown in Fig. 14. 

We confirm that a PLE (or PDE) with a constant evo- 
lution parameter is a good approximation for the redshift 
range z G [0.7,1.7]. At larger redshift, the results of our 
analysis of the LBQS are consistent with a maximum lu- 
minosity (or density) around z = 2.5, but the reliability 
of these results is not yet established due to likely selec- 
tion biases. The determination of this k(z) dependence 
towards larger z is essential for the understanding of the 
birth and growth of quasars at high redshifts, and the re- 
lation of the quasar phenomenon with star bursts in the 
primordial universe. 

Appendix 

A rough estimation of the errors on k^ and kjj has been 
made, using the correlation we found between (V/V max ) 
and fci or ku values. (V/V max } has been calculated for 
each bin making the hypothesis of a zero evolution, k^ ver- 
sus (V/Vmax) is shown in Fig. f4. The dashed line and the 
squares correspond to binning (Bf ) of AAT catalogue. The 
dashed-dotted line and the stars correspond to AAT's bin- 
ning (B2). The solid line and the triangles are the LBQS 
results. On the upper left we show a typical error bar on 
the (V/V max ) estimated as <J(v/v m 

The coefficients of the linear approximation fc_L pm , = 
at + [3l (V/V max ) are given in Table 3 (also for ko), for 
AAT (binnings (BI) and (B2)) and LBQS. We notice that 
there a unique correlation may be adopted between fei , 
and (V/V max ), as the error bars are much more important 
than the variations of the coefficients and {3l- Fig. f4 
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Fig. 14. kh versus (V/Vmax) in the case of a power law 
parametrization. , □: AAT's first binning. — , *: AAT's sec- 
ond binning. — , A: LBQS. All calculations have been made for 
the cosmological model (fio, A) = (0.5,0.5). In the upper-left 
corner, a typical error bar CT(v/v ma x) = 1/V12JV. 



clearly illustrates that if the quasar population docs not 
evolve (fci or fcjj = 0), then their (V/V max ) is equal to 
0.5, as expected. 



Table 3. Numerical values for the coefficients in the 
linear approximation fei, , = «i + (3l {V/Vmax} and 
kD pw i = cto + (3d (V/Vmax) in a PLE and a PDE hypothe- 
sis, respectively. 



Catalogue/ 


CtL 


(3l 


an 


(3d 


binning 










A AT (Bl) 


-35 


70.29 


-22.95 


47.75 


AAT (B2) 


-41.97 


83.12 


-36.74 


73.08 


LBQS 


-36.6 


72.26 


-51.67 


102.59 
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In the case of an exponential parametrization and in 
order to define the linear correlation, we must take under 
consideration the differentiation of the evolution parame- 
ters with cosmological models. 

The existence of a correlation between these two quan- 
tities was somehow expected, in the sense that the higher 
the value of (V/Vmax) is under a zero evolution hypoth- 
esis, the higher k L (or fcjj) must be in order to insure a 
(V/V max ) = 0.5 in a PLE (or PDE) hypothesis. 
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